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rObservational work conducted over the last few decades indicates that 
all massive galaxies have supermassive black holes at their centres. 
F^^lthough the luminosities and brightness fluctuations of quasars in 
ftle early Universe suggest that some are powered by black holes with 
glasses greater than 10 billion solar masses ' ^, the remnants of these 
M>bjects have not been found in the nearby Universe. The giant ellip- 
^^cal galaxy Messier 87 hosts the hitherto most massive known black 
Cllole, which has a mass of 6.3 billion solar masses ^ Here we report 
'""tliat NGC 3842, the brightest galaxy in a cluster at a distance from 
Ij ^arth of 98 megaparsecs, has a central black hole with a mass of 9.7 
pillion solar masses, and that a black hole of comparable or greater 
Q^ass is present in NGC 4889, the brightest galaxy in the Coma cluster 
|~-^at a distance of 103 megaparsecs). These two black holes are sig- 
f^ificantlv more massive than predicted by linearly extrapolating the 
T-^idely-used correlations between black hole mass and the stellar ve- 
^vlpcity dispersion or bulge luminosity of the host galaxy ' Although 
tftf'Kf correlations remain useful for predicting black hole masses in 
^— l^ss massive elliptical galaxies, our measurements suggest that differ- 
^— ent evolutionary processes influence the growth of the largest galaxies 
Hand their black holes. 

? Empirical scaling relations between black hole mass (Mbh), galaxy 
fr^^ulge velocity dispersion (a), and luminosity (L) are commonly used to es- 
^mate black hole masses, because for most galaxies we are unable to make 
a direct measurement. Estimates of the number density of black holes in a 
given mass range thus depend upon the empirically determined A/bh — cr 
and Mbh — L relations over an appropriate range of galaxy masses. Di- 
rectly measuring A/bh from the kinematics of stars or gas in the vicinity of 
the black hole is particularly difficult at the highest galaxy masses, because 
massive galaxies are rare, their typical distances from Earth are large, and 
their central stellar densities are relatively low. The most massive galaxies 
are typically brightest cluster galaxies (BCGs), that is, giant ellipticals that 
reside near the centres of galaxy clusters. 

We have obtained high-resolution, two-dimensional data of the line- 
of-sight stellar velocities in the central regions of NGC 3842 and NGC 
4889 using integral field spectrographs at the Gemini North and Keck 2 
telescopes, in Hawaii. The stellar luminosity distribution of each galaxy 
is provided by surface photo metry from NASA's Hubble Space Telescope 
and ground-based telescopesCM NGC 3842 is the BCG of Abell 1367, a 
moderately rich galaxy cluster. NGC 4889 is the BCG of the Coma cluster 
(Abell 1656), one of the richest nearby galaxy clusters. We targeted these 



two galaxies because they have relatively high central surface brightnesses 
and lie at an accessible distance for direct measurements of A/bh. 

We measured the distribution of stellar velocities at 82 different loca- 
tions in NGC 3842. The line-of-sight velocity dispersion in NGC 3842 
is between 270 and 300 km s^^ at large galactocentric radii (r) and rises 
in the central 0.7 arcsec (r < 330 pc), peaking at 326 km s^^ (Figs [T| 
[2](. We determined the mass of the central black hole by constructing 
a series of orbit superposition models Each model assumes a black 
hole mass, stellar mass-to-light ratio (A/t/L) and dark matter profile, and 
generates a library of time-averaged stellar orbits in the resulting gravi- 
tational potential. The model then fits a weighted combination of orbital 
line-of-sight velocities to the set of measured stellar velocity distributions. 
The goodness-of-fit statistic is computed as a function of the assumed 
values of A/bh and the stellar mass-to-light ratio. Using our best-fitting 
model dark matter halo, we measure a black hole mass of 9.7 X 10^ solar 
masses (A/0), with a 68% confidence interval of 7.2 - 12.7 x 10^ Mq. 
Models with no black hole are ruled out at the 99.996% confidence level 
( Ax^ = 17. 1). We find the stellar mass-to-light ratio to equal 5. 1 A/q /Lq 
in R band (Lq, solar luminosity), with a 68% confidence interval of 
4.4A/0/L0-5.8A/Q/L0). 

We measured stellar velocity distributions at 63 locations in NGC 4889 
and combined our measurements with published long-slit kinematics at 
larger radii The largest velocity dispersions in NGC 4889 are located 
across an extended region on the east side of the galaxy. The stellar orbits 
in our models are defined to be symmetric about the galaxy centre, so we 
constrain A/bh by running separate trials with velocity profiles from four 
quadrants of the galaxy. The best-fitting black hole masses from the four 
quadrants range from 9.8 x 10^ A/© to 2.7 x 10^° Mq. All quadrants favor 
tangential orbits near the galaxy centre, which cause the line-of-sight ve- 
locity dispersion to decrease even as the internal three-dimensional veloc- 
ity dispersion increases toward the black hole. Although no single model 
is consistent with all of the observed kinematic features in NGC 4889, we 
can define a confidence interval for A/bh by considering the most extreme 
confidence limits from the cumulative set of models. The corresponding 
68% confidence interval is (0.6 - 3.7) x 10^° Mq- We adopt a black hole 
mass of 2.1 x 10^° Mq, corresponding to the midpoint of this interval. 

Figure [3] shows the A/bh — o" and A/bh — L relations, using data 
compiled from studies published before the end of August 201 1, plus our 
measurements of NGC 3842 and NGC 4889. Tabulated data with refer- 
ences are provided as Supplementary Information. The most widely used 
form for both relations is a power law with a constant exponent. Straight 
lines in Fig.[3]show our fits to A/bh(o") and Mbh{L). The relationship 
between a and L, however, flattens at high galaxy masses, and constant- 
exponent power laws for the A/bh — cr and A/bh — L relations produce 
contradictory predictions for A/bh in this mass range Direct measure- 
ments of A/bh in higher mass galaxies will compel the revision of one or 
both of the A/bh — rr and A/bh — L relations. 
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The average velocity dispersion in NGC 3842 is 270kms~^, mea- 
sured outside the black hole radius of influence (1.2 arcsec or 570 pc) and 
inside the two-dimensional half-light radius (38 arcsec or 18 kpc). Al- 
though NGC 3842 hosts a black hole more massive than any previously 
detected, its average dispersion ranks only fourteenth among 65 galaxies 
with direct measurements of A/bh- Its luminosity ranks fifth in this sam- 
ple of galaxies and is exceeded only by other BCGs. On the basis of a 
and L for NGC 3842, our revised Mbh — cr and Mbh — L relations pre- 
dict A^BH = 9.1 X 10* A'/© and 2.5 x 10^ Mq, respectively. Similarly, for 
NGC 4889 the respective predictions are 3.3 x 10^ Mq and 4.5 x lO'' Mq- 
These predictions are smaller than our direct measurements of A^bh, by 
1.6-4.6 times the 1-s.d. scatter in the A^bh — o" and A/bh — L relations 
Four measurements of A/bh in BCGs existed before this work. Two 
measurements based on gas dynamics and one based on stellar dynamics 
all lie within 1.2 s.d. of our revised fits to the A/bh — o" and A/bh — L 
relations E32i. Yet the measurement of A/bh in NGC 1316, the BCG of 
the Fornax cluster, is 3.4 s.d. less than that predicted by our A/bh — L 
relation The high scatter indicated by this collection of measurements 
reveals large uncertainties in the standard practice of using galacitc a or L 
as a proxy for the central black hole mass in giant elliptical galaxies and 
their predecessors. 

Several BCGs within 200 Mpc of Earth are at least twice as lumi- 
nous as NGC 3842, and three times as luminous as M87, which hosted 
the most massive black hole known before this work. In spite of their ex- 
treme luminosities, BCGs have velocity dispersions similar to those of the 
most massive field elliptical galaxies. Yet the most massive black holes are 
found predominantly in BCGs (Fig. [3]l. How galaxies are assembled and 
the role of gas dissipation affect the correlations (or lack thereof) among 
A/bh, o", and L. Simulations of mergers of gas-rich disk galaxies are able 
to produce remnant galaxies that follow the observed A/bh — o" correla- 
tion in Fig. 3a, over the intermediate mass range A/bh « 10^ - 10^ A/q 
(refs 18, 19). By contrast, simulated mergers of elliptical galaxies with 
low-angular momentum progenitor orbits increase A/bh and L by similar 
numerical factors, without increasing the velocity dispersionl^^ Because 
these mergers are a likely path to forming the most massive galaxies, the 
A^bh — correlation may steepen or disappear altogether at the highest 
galaxy masses. Massive elliptical galaxies retain residual quantities of gas 
even after the decline of star formation. Accretion of this gas onto the 
galaxies' central black holes could help increase A/bh and further steepen 
the Mbh — cr and A^bh — L relations. 

Black holes in excess of 10^° A/q are observed as quasars in the early 
Universe, from 1.4 x 10^ to 3.3 x 10® yr after the Big Bang ^ (redshift, 
2 = 2 — 4.5). Throughout the last 1.0 x 10^° yr, however, these extremely 
massive black holes have not been accreting appreciably, and the average 
mass of the black holes powering quasars has decreased steadily. Quasar 
activity and elliptical galaxy formation are predicted to arise from similar 
merger-triggered processes, and there is growing evidence that present-day 
massive elliptical galaxies once hosted the most luminous high-redshift 
quasars Yet definitive classification of these quasars' host galaxies has 
remained elusive. 

Our measurements of black holes with masses of around 10^° A/q 
in NGC 3842 and NGC 4889 provide circumstantial evidence that BCGs 
host the remnants of extremely luminous quasars. The number density of 
nearby BCGs (~ 5 x 10^^ Mpc^^) is consistent with the number den- 
sity of black holes (~ 3 x lO"'^ to lO"'' Mpc"^) with masses between 
10'' Mq and 10^° Mq predicted from the A/bh - L relation and the lu- 
minosity function of nearby galaxies. Furthermore, both quantities agree 
with predictions based on the black hole masses and duty cycles of quasars. 
The black hole number density predicted from the A/bh — o" relation, how- 
ever, is an order of magnitude less than the inferred quasar 

populationlHISl. 

These two predictions can be reconciled if the Mbh — cr relation has up- 
ward curvature or a large degree of intrinsic scatter in A/bh at the high- 
mass end, as suggested by our new measurements. With improvements in 



adaptive optics instrumentation on large optical telescopes and very-long 
baseline interferometry at radio wavelengths, black holes are being sought 
and detected in increasingly exotic host galaxies. Along with our measure- 
ments of the black hole masses in NGC 3842 and NGC 4889, future mea- 
surements in other massive galaxies will quantify the cumulative growth 
of supermassive black holes in the Universe's densest environments. 
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Figure 1 Two-dimensional maps of the line-of-sight stellar velocity dispersions. 

The maps show the central regions of NGC 3842 (a) and NGC 4889 (b) ob- 
served using the GMOS spectrograph on the 8-m Gemini North telescope. 
Additional l^inematics at large radii were measured using the VIRUS-P spectro- 
graph ^'^ at the 2.7-m Harlan J. Smith telescope, and additional high-resolution 
data were acquired with the OSIRIS spectrograph at the 10-m Keck 2 tele- 
scope. GMOS, OSIRIS, and VIRUS-P are all Integral field spectrographs, 
which record spectra at multiple positions In a two-dimensional spatial array. 
The horizontal dashed line in each panel traces the major axis of the galaxy. 
The median 1 s.d. errors in velocity dispersion are 12 km s^^ and 20 km s^^ 
for NGC 3842 and NGC 4889, respectively. In NGC 4889 the highest velocity 
dispersions, near 410 km s^^, are located on the east side of the galaxy, at 
least 1.1 arcsec from the centre. 
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Figure 2 One-dimensional profiles of the line-of-sight velocity dispersions, (a) 

Dispersion versus radius in NGC 3842, after averaging data at a given ra- 
dius, based on measurements with GMOS (black circles) and VIRUS-P (red 
diamonds). The solid blue line Is the projected line-of-slght dispersion from 
our best-fitting stellar orbit model of NGC 3842. (b) Dispersion versus radius 
along the major axis of NGC 4889, measured from GMOS (black circles) and 
the William Herschel Telescope (green diamonds). The maximum velocity 
dispersion occurs at r = 1.4 arcsec. The solid blue line is the projected line- 
of-sight dispersion from our best-fitting orbit model using data from the east 
side of NGC 4889 (r > 0). The dashed orange line is from our best-fitting 
orbit model using data from the west side of NGC 4889 {r < 0). Error bars, 1 
s.d. 
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Figure 3 Correlations of dynamically measured black hole masses and bulk properties of host galaxies, (a) Black hole mass, Mbh, versus Stellar velocity 
dispersion, ct, for 65 galaxies with direct dynamical measurements of Mbh- For galaxies with spatially resolved stellar kinematics, a is the luminosity-weighted 
average within one effective radius (Supplementary Information), (b) Black hole mass versus V-band bulge luminosity, Lv {Lq^v, solar value), for 36 early- 
type galaxies with direct dynamical measurements of Mbh- Our sample of 65 galaxies consists of 32 measurements from a 2009 compilation ^, 16 galaxies 
with masses updated since 2009, 15 new galaxies with Mbh measurements, and the two galaxies reported here. A complete list of the galaxies is given in 
Supplementary Table 4. BGGs (defined here as the most luminous galaxy in a cluster) are plotted in green, other elliptical and SO galaxies are plotted in red, 
and late-type spiral galaxies are plotted in blue. The black hole masses are measured using dynamics of masers (triangles), stars (stars), or gas (circles). 
Error bars, 68% confidence intervals. For most of the maser galaxies, the error bars in Mbh are smaller than the plotted symbol. The solid black line in (a) 
shows the best-fitting power law for the entire sample: log (Mbh/ Mq) = 8.29-1-5.12 log^(,[<T/(200kms^^)]. When early-type and late-type galaxies 
are fit separately, the resulting power laws are logio(AfBH/ Mq) = 8.38-1- 4.53 logio[o-/(200kms^^)] for elliptical and SO galaxies (dashed red line), and 
logio(MBH/ Mq) = 7.97 -I- 4.581ogio[CT/(200kms^^)] for spiral galaxies (dotted blue line). The solid black line in (b) shows the best-fitting power law: 
logio(MBH/ Mq) = 9.16-1- 1.16 logio(Lv/10^^Lo)- We do not label Messier 87 as a BOG, as is commonly done, as NGC 4472 in the Virgo cluster is 
0.2 mag brighter 
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Supplementary Information 

In the first section, we describe spectroscopic data of NGC 3842 and 
NGC 4889 and our procedures for measuring stellar kinematics. In the 
second section, we describe photometric data. In the third section we com- 
pare overlapping kinematic measurements from GMOS, OSIRIS, VIRUS- 
P, and long-slit spectroscopy. The fourth section contains a summary of our 
stellar orbit models. In the fifth section we describe the unusual kinematics 
and stellar orbital structures in NGC 4889, and justify our measurement of 
Mbh- In the final section we describe our updated sample and fits to the 
Mbh — cr and Mbu — L relations. 

We tabulate our spectroscopic observations, our results from modeling 
each quadrant of NGC 4889, our fits to the Mbh — o" and Mbh — L rela- 
tions, and the galaxy sample used to fit the relations. Our measured line- 
of-sight velocity distributions (LOSVDs) for NGC 3842 and NGC 4889 
are available as online data. 

Spectroscopic Data and LOSVD Extraction 

We map stellar orbital motions in NGC 3842 and NGC 4889 by measuring 
the LOSVD for different regions in each galaxy. Each LOSVD is deter- 
mined by fitting a composite template stellar spectrum to a fully reduced 
spectrum of the galaxy. The LOSVDs are non-parametric probability dis- 
tributions, defined at each of 15 velocity bins. We use a Maximum Penal- 
ized Likelihood technique to optimize the LOSVD value in each velocity 
bin while si multaneously optimizing the weights of individual template 

starsinniEzi. 

In Table[T] we summarize our observations with the integral field spec- 
trographs GMOS, OSIRIS, and VIRUS-R Figures |4|6] illustrate a sample 
galaxy and template spectrum from each instrument. Our data from GMOS 
only cover radii within 3.8 arcseconds (1.8 kpc) of the center of NGC 3842 
and NGC 4889, and by themselves cannot fully remove the degeneracies 
between Mbh and Af*/L. For NGC 3842, our VIRUS-P measurements 
cover radii out to 35.3 arcseconds (16.8 kpc) and can distinguish the en- 
closed stellar mass profile from the galaxy's dark matter halo. This allows 
for an accurate determination of A/*/L, such that the GMOS data can ac- 
curately constrain Mbh- At radii from 3.6 to 23.0 arcseconds (1.8 — 11.5 
kpc) along the major axis of NGC 4889, we use Gaussian velocity profiles 
from Loubser et al. (2008)^^. 

Our GMOS spectra for NGC 3842 and NGC 4889 are centered on the 
calcium triplet absorption lines near 860 nm. A sample GMOS spectrum 
for each galaxy is shown in Figure |4] demonstrating the clean line profiles 
that are typical for this spectral region. Another advantage to using the 
calcium triplet is that kinematic measurements are not highly sensitive to 
the stellar template usedEU. 

OSIRIS data of NGC 3842 were acquired with the 0.05 arcsecond 
spatial scale and the broad i/-band filter, which spans a large number of 
atomic and molecular absorption features. To measure kinematics, we fit 
carbon monoxide band heads at 1598, 1619, 1640, and 1661 nm, and a 
deep magnesium feature near 1500 nm. The most severe source of noise 
in our OSIRIS spectra is residual narrow-line emission from the night sky. 
This background emission varies rapidly and is only partially corrected by 
recording 15 -minute sky frames in between pairs of 15-minute science ex- 
posures. Even when masked, the contaminated channels represent a non- 
negligible loss of spectral information. This loss must be countered by 
increasing the signal-to-noise ratio in the usable parts of the spectrum. To 
achieve adequate signal-to-noise, we bin the data to the same spatial re- 
gions as the overlapping LOSVDs from GMOS. Each of the final bins 
contains approximately 80 OSIRIS spatial pixels. 

VIRUS-P data were acquired in low-resolution mode, which provides 
broad wavelength coverage. The poorest instrumental resolution over our 
field is 0.56 nm full width at half-maximum (FWHM), which is adequate 
considering that stellar velocity dispersions in NGC 3842 are typically 
near 250 kms^^, corresponding to a resolution of approximately 0.7 nm 
FWHM x(A / 360 nm). We fit VIRUS-P spectra with 16 template stars 



from the Indo-US library'^, spanning spectral types from B9 to M3 and in- 
cluding stars with sub-solar and super-solar metallicities. However, some 
spectral regions require additional adjustments to account for metallicities 
and elemental abundance ratios outside the range of our template library. 
We do not attempt to fit the entire spectral range simultaneously, but in- 
stead follow the same fitting procedure used by Murphy et al. (201 1) for 
VIRUS-P data of M87. Each spectrum is divided into four sub-regions, as 
illustrated in Figure |6] We independently determine the best-fit LOSVD 
for each sub-region and discard any sub-regions that fail to produce a be- 
lievable fit. We then average the LOSVDs derived from individual spectral 
sub-regions. 

Our stellar orbit models are axisymmetric, and so each of our final 
LOSVDs must represent an average over four quadrants of the galaxy. In 
order to preserve any rotational signal along the major axis, we invert the 
velocities from LOSVDs on the south side of NGC 3842 and the west 
side of NGC 4889. However, neither galaxy shows strong rotation. In 
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Figure 4 | GMOS spectra of NGC 3842 and NGC 4889. Each spectrum cor- 
responds to the center of the galaxy (r < 0.25 arcseconds). The upper 
spectrum is NGC 3842 (black), overlaid with the best-fitting, LOSVD-convolved 
template spectrum (thick, light blue). The middle spectrum is NGC 4889, over- 
laid with the best-fitting template spectrum. The dotted portion of the spectrum 
was excluded from the LOSVD fitting. The lower spectrum is template star HD 
73710 (G9III), before convolution with the LOSVD. 



X 1.1 
s 

E 1.0 

I "-^ 

5 0.8 

B 0.7 

9 

I 0.5 





1500 



1550 1600 1650 



Figure 5 | OSIRIS spectrum of the central region of NGC 3842 (r < 0.25 arc- 
seconds, or 1 10 pc). The upper spectrum is NGC 3842 (black), overlaid with 
the best-fitting, LOSVD-convolved template spectrum (thick, light blue). The 
red dashed lines in the galaxy spectrum are residuals from imperfectly sub- 
tracted sky emission. The lower spectrum is the best-fit composite template 
before convolution with the LOSVD. Our observed template stars fit spectra of 
NGC 3842 poorly across the dotted region from 1510 to 1590 nm, and there- 
fore this region is excluded from kinematic fitting. 



NGC 3842, the resulting kinematics are sufficiently symmetric to average 
the LOSVDs from opposite sides of the galaxy. In NGC 4889, we have 
modeled four quadrants of the galaxy independently. 

Our stellar orbit models require an estimate of the point spread func- 
tion (PSF) for each instrument at the time of spectroscopic observations. 
We estimate the PSF for our GMOS data for both NGC 3842 and NGC 
4889 from wide-field images taken during target acquisition. Although 
the models discussed herein assume a 0.4-arcsecond PSF for GMOS data, 
we have run a small number of models with a 0.7-arcseond PSF and have 
found no significant changes in our results. While observing NGC 3842, 
we switched from the OSIRIS spectrograph to the OSIRIS imaging cam- 
era every few hours, and observed the adaptive optics tip/tilt star. We can 
tolerate a large degree of uncertainty in our measured PSF, as we have re- 
binned OSIRIS data of NGC 3842 to spatial scales that are several times 
coarser than the diffraction-limited FWHM. 

Photometric Data 

Our stellar orbit models are constrained to reproduce the observed stel- 
lar light profile of each galaxy, which requires accurate measurements of 
each galaxy's surface brightness profile over a large radial range. For radii 
out to 10 arcseconds, we adopt high-resolution J-band (800 nm) surface 
brightness profiles, obtained with WFPC2 on the Hubble Space Telescope, 
and deconvolved with the instrumental PSF At larger radii out to 115 
arcseconds, we use i?-band (600 nm) data obtained with the 2.1 m tele- 
scope at Kitt Peak National Observatory (KPNO). The KPNO data have a 
field of view of 5.2 x 5.2 arcminutes, which enables accurate sky subtrac- 
tion. We combine the individual profiles from WFPC2 and KPNO data at 
overlapping radii between 5 and 10 arcseconds, accounting for the average 
R — I color over these radii. To compute the luminosity density profile 
of each galaxy, we deproject the surface brightness profile while assuming 
spheroidal isodensity contours'^. 



Kinematic Data from Different Instruments 

Stellar kinematics along the major axes of NGC 3842 and NGC 4889 were 
previously measured by Loubser et al. (2008) using the ISIS long- 
slit spectrograph on the William Herschel Telescope (WHT). The WHT 
measurements of NGC 3842 are in good agreement with our GMOS mea- 
surements for r < 3.0 arcseconds. The WHT data for NGC 4889 agree 
with our GMOS measurements on the west side of NGC 4889 and at radii 
between 1 and 2.4 arcseconds on the east side. However, they do not re- 
produce our measurement of the large central drop in stellar velocity dis- 
persion (Figure |7j. For a more direct comparison, we have rebinned our 
integral-field spectra to match the 0.4-arcsecond spatial sampling and 1 .0- 
arcsecond slit width of the WHT data. We find that rebinning alleviates the 
velocity dispersion discrepancy for all but the central point (r = 0). Our 
detection of a significantly sharper decrease in velocity dispersion is con- 
sistent with the superior seeing conditions of our data (0.4 arcseconds for 
GMOS, versus 1.0 arcseconds for WHT). Overall, we find our kinematic 
measurements in NGC 3842 and NGC 4889 to be broadly consistent with 
the independent measurements by Loubser et al., indicating that our kine- 
matic extraction method has low systematic errors. 

For N3842, our VIRUS-P measurements are also consistent with WHT 
measurements, which extend along the major axis to r = 20.8 arcseconds. 
We prefer using data from VIRUS-P in our stellar orbit models because 
they extend to larger radii and provide full two-dimensional spatial sam- 
pling. 

Near the center of NGC 3842, OSIRIS and GMOS provide inde- 
pendent measurements of stellar kinematics. We have binned data from 
OSIRIS and GMOS at identical spatial scales out to r = 0.7 arcseconds 
(330 pc) and have run orbit models fitting LOSVDs from OSIRIS and 
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Figure 6 | VIRUS-P spectrum of NGC 3842, corresponding to a semi-annulus 
with an inner radius of 17.0 arcseconds (7.7 kpc) and an outer radius of 24.5 
arcseconds (1 1 .0 kpc). The four panels contain different sub-regions of the 
galaxy spectrum. Each sub-region is evaluated independently for a best- 
fit LOSVD and best-fit composite template spectrum. The upper spectrum 
in each panel is NGC 3842 (black), overlaid with the best-fitting, LOSVD- 
convolved template spectrum (thick, light blue). Dotted portions of the galaxy 
spectrum have been masked from the fit. The lower spectrum in each panel is 
the best-fit composite template before convolution with the LOSVD. 
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Figure 7 | (a) Radial velocity vs. radius in NGC 4889. (b) Velocity dispersion 
vs. radius. Green triangles are measurements from Loubser et al. (2008)1^, 
using the ISIS spectrograph on the William Herschel Telescope (WHT). Black 
circles are our measurements using the GMOS intrgral-field unit. Error bars 
represent one standard deviation. On the west side and at radii from 1.0 to 
2.4 arcseconds on the east side, the GMOS and WHT measurements agree 
within errors. At r < 1.0 arcseconds on the east side, only GMOS detects a 
local minimum in velocity dispersion. The central feature might be unresolved 
in the WHT data due to worse seeing (1.0 arcseconds). 



GMOS simultaneously (as well as VIRUS-P data at large radii). Including 
the OSIRIS data causes the best-fit value of Mbh to decrease by up to 
23%, and the best-fit value of h'h/L to increase by as much as 8%. This 
occurs because OSIRIS data show a less drastic increase in velocity dis- 
persion than data from GMOS. In spite of these differences, results with 
and without OSIRIS data are consistent at the 68% confidence level. 

Models fitting OSIRIS and GMOS data together yield higher aver- 
age values per LOSVD. This is true even if we ignore the central re- 
gions where LOSVDs from OSIRIS and GMOS are not fully consistent. 
Even with several template stars, the overlapping absorption features in 
the H-band spectral region are difficult to model, and the LOSVDs derived 
from OSIRIS data may have systematic errors. Consequently, we judge 
the models with only GMOS and VIRUS-P data to be more reliable. 

Stellar Orbit Models and Statistical Analysis 

We generate models of NGC 3842 and NGC 4889 using Schwarzschild's 
method^, in which test particle orbits are computed in a static axisym- 
metric gravitational potential. We assume that each galaxy contains three 
mass components: stars, a central black hole, and an extended dark matter 
halo. The stellar mass density is assumed to follow the same profile as 
the observed luminosity density, with a constant stellar mass-to-light ratio, 
Mi,/L. Our modeling procedures are similar to those for NGC 6086, de- 
scribed in McConnell et al. (201 1) 

Each orbit in the model is assigned a scalar weight, and the set of best- 
fit orbital weights is determined by comparing projected LOSVDs from 
the orbits to the observed LOSVDs for the galaxy. Each observed LOSVD 
spatially maps to a linear combination of bins within the model, according 
to the spatial boundaries of the corresponding spectrum and the PSF of 
the observations. A corresponding model LOSVD is computed from the 
projected velocity distributions of individual orbits in each spatial bin, the 
appropriate combination of spatial bins, and the orbital weights. The best- 
fit weights are determined by the method of maximum entropy with the 
fixed constraint that the summed spatial distribution of all weighted orbits 
must match the observed luminosity density profile. The essential output 
of each model is a measurement of x^, which defines the goodness of fit 
between our observed LOSVDs and the model LOSVDs, using the opti- 
mal combination of orbital weights. We determine the best-fit values and 
confidence intervals in A/bh and K'h/ L by evaluating the relative likeli- 
hood between models with different assumed values of Mbh and M-^jL. 
Figure [s] illustrates the behavior of with respect to Mbh and M^,/L, 
for models with our best fitting dark matter halo. 

Models of NGC 4889 

Stellar kinematics in NGC 4889 are asymmetric with respect to the major 
and minor axes of the galaxy. Integral-field data from GMOS reveals ve- 
locity dispersions above 410kms~^ on the east side of the galaxy, while 
the velocity dispersion rarely exceeds 380kms~^ on the west side. This 
asymmetry prevents NGC 4889 from being fully described by a single set 
of axisymmetric orbit models. In order to place upper and lower bounds 
on the central black hole mass, we have run four suites of models, each fit- 
ting kinematics from one projected quadrant of NGC 4889. The northeast, 
southeast, and northwest quadrants yield consistent black hole masses, 
spanning a 68% confidence interval of A^bh = 1-0 - 3.7 x 10^° Af©. 
The southwest quadrant has a maximum velocity dispersion of 373 km/s 
and yields a 68% confidence interval of Mbh = 0.6 - 1.7 x 10^" Mq. 

We have run an additional set of models to approximate A/bh in the 
case of an off-center black hole. We apply a constant spatial offset of 1.4 
arcseconds (700 pc) to the kinematics on the east side of the galaxy, such 
that the highest velocity dispersion is aligned with the center of the model 
gravitational potential. These models cannot be fully trusted because the 
kinematic and photometric data are misaligned. Still, the resulting 68% 
confidence interval for A/bh falls entirely within the range bracketed by 
the models from different quadrants. Results from individual trials are 



listed in Table [2] Figure [9] illustrates versus A^bh for each series of 
models, after marginalizing over Mi,/ L. 

Although three-dimensional stellar velocities must increase in the 
vicinity of a black hole, a deficiency of radial orbits can produce a cen- 
tral minimum in the line-of-sight velocity dispersion, as we observe in 
NGC 4889. Indeed, our best-fitting models of NGC 4889 exhibit tangen- 
tial bias at small radii, as shown in Figure[To] In contrast, models without 
a black hole reproduce the central drop in velocity dispersion with a nearly 
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Figure 8 | Map of versus A/bh and AU/L for stellar orbit models of NGC 
3842. The models fit data from GMOS and VIRUS-P. The diagonal trend in 
Indicates the degeneracy between stellar mass and black hole mass near 
the center of NGC 3842. For two free parameters with Gaussian likelihood 
distributions, the 68% confidence interval Is defined where x^-xLin < 2.30, 
illustrated by the thick white contour. We obtain 68% confidence intervals of 
7.2 - 12.7 X W'^ Mq for A//bh, and 4.4 - S.SMq/Lq for M^,/L. The 
median values, which we adopt as our final measurements, are A/bh = 9-7 x 
10^ Mq and M^,/L = 5.1 Mq/Lq. 




Figure 9 | x vs. A/bh for NGC 4889, after marginalizing over M^/L. Each 
line with symbols represents a models constrained by different set of LOSVDs 
from GMOS. Black squares, purple circles, red diamonds, and green triangles 
each use LOSVDs from a different quadrant of the galaxy. Blue stars represent 
models with spatially offset LOSVDs, to match the largest velocity dispersion 
with the center of the gravitational potential. Considering all models, the 68% 
confidence Interval for A/bh is 0.6 - 3.7 x 10^° Mq. 
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isotropic orbital distribution. However, these models yield a worse overall 
fit, indicated by higher values of in each quadrant. Models of NGC 
3842 exhibit a similar but less severe trend, consistent with the modest 
increase in line-of-sight velocity dispersion toward the center. 

The individual quadrants of NGC 4889 represent large variations in 
stellar kinematics, but each quadrant still partially constrains the enclosed 
mass within the central few arcseconds. By adopting the most extreme 
range of confidence limits, Mbh ~ 0.6 — 3.7 x 10^" Mq, we only ex- 
clude black holes whose gravitational influence would contradict our entire 
field of data. Further extensions to this confidence interval should only re- 
flect overall systematic biases. Large systematic biases in our kinematic 
measurements are unlikely, as demonstrated by their agreement with inde- 
pendent measurements by Loubser et al. (2008) Our models of NGC 
4889 assume an edge-on inclination. This is indirectly supported by the 
observed axis ratio of 0.7, which implies a relatively eccentric intrinsic 
shape even for an edge-on system. Models with a more face-on inclination 
might yield a systematically higher black hole mass A fundamental 
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Figure 10 | Orbital anisotropy in models of NGC 3842 and NGC 4889. Three- 
dimensional stellar velocities in the models are divided into radial and tan- 
gential components. The ratio of velocity dispersions, o-rad/o"tan, varies with 
radius and the assumed Mbh in the model, (a) The thick solid line repre- 
sents the best-fitting model for NGC 3842, with A/bh = 8.5 x 10^ Mq. 
The thin dashed line represents the best model with no black hole, (b) Thick 
solid lines represent the best-fitting model for each quadrant in NGC 4889, 
with A/bh = 1.3 X 10^"Afo (northeast), 2.8 x 10^° A/q (northwest), 
2.4 X 10^° Mq (southeast), and 6.5 x 10*' Mq (southwest). Thin dashed 
lines represent models with no black hole. For NGC 3842 and each quadrant 
of NGC 4889, our kinematic data are best fit with a massive black hole and a 
prevalence of tangential orbits at r < 2 arcseconds (1 kpc). In NGC 4889, 
this tangential bias is responsible for the observed decrease in line-of-sight 
velocity dispersion near the center. 



assumption of all orbit superposition models is that the stellar motions re- 
flect a steady-state gravitational potential, rather than transient conditions. 
These models could misrepresent the range of allowed black hole masses if 
the observed kinematics in NGC 4889 reflected a temporary phenomenon 
such as an ongoing galaxy merger. NGC 4889 appears photometrically 
undisturbed, reducing the likelihood of such an event. 

Because our adopted confidence interval places large error bars on 
Mbh in NGC 4889, this galaxy has relatively little weight in our fits to 
the A/bh — o" and A/bh — i relationships. A systematic error in our mea- 
surement would produce a minimal bias in the best-fit relations. Likewise, 
our discussion of steepening and scatter at the high-mass ends of the cor- 
relations depends upon several objects and is not highly sensitive to the 
measurement in NGC 4889. 



Power-Law Fits to the Mbh - o and Mbh - 1^ Relations 

We revisit the A/bh — o" and Mbh — i relations by updating the sample 
of 49 black holes from Gultekin et al. (2009; hereafter G09) which was 
compiled from earlier studies. Including NGC 3842 and NGC 4889, we 
add 17 galaxies with recently measured black hole masses to the sample. 
The new objects include two more BCGs ei ght a ctive galactic nu- 
clei with high-precision maser-based measurements"^^] and two galaxies 
with pseudobulges We also include updated black hole masses for 16 
other galaxies in the 2009 sample, based on ste llar orbi t models with dark 
matter halos and more thorough orbit libraries 0^3^311. in particular, the 
revised masses for M60 and M87 are twice as large as the earlier values. 

Our updated sample uses the same selection criteria as G09; in partic- 
ular, only direct dynamical measurements of Mbh are included. G09 esti- 
mated galaxy distances by assuming a Hubble parameter i/o ~ 70 km s^^ 
Mpc^^, and rescaled their sample of black hole masses accordingly 
(A/bh oc D). We have followed this convention for NGC 3842, NGC 
4889, and the rest of our sample. Our fits to A/bh (o") do not include upper 
limits. Updated models of one galaxy, NGC 2778, do not produce a signif- 
icant black hole detection; after removing this object, our updated sample 
contains 65 black hole masses. 

We define a in the same manner as the 2009 sample. Wherever pos- 
sible, we use the luminosity-weighted effective velocity dispersion, mea- 
sured using spatially resolved data out to one effective radius. G09 found 
no evidence of systematic bias between this definition of a and more ubiq- 
uitous single-aperture measurements. Nonetheless, several galaxies in G09 
have measurements of a that include data at very small radii, within which 
the central black hole directly influences the stellar velocity dispersion. 
This is inappropriate for studies that wish to treat A/bh and a as funda- 
mentally independent variables. We have therefore re-evaluated a in three 
galaxies with large black hole masses and available spatially resolved kine- 
matics, by excluding data within the black hole's radius of influence. We 
use this same treatment to measure a for NGC 3842 and NGC 4889. For 
M87, we adopt the updated values of A/bh and a from Gebhardt et al. 
(2011) Following G09, we assume all measurements of a have an un- 
certainty of at least 5%. 

To fit the A/bh — L relationship, we only consider early-type galaxies 
for which the stellar luminosity of the spheroidal component can be mea- 
sured reliably. Including NGC 3842, we add 6 galaxies to the A/bh — L 
sample of G09. We exclude NGC 2778. We also exclude NGC 3607 
and NGC 4564, for which the literature contains large discrepancies in the 
measured luminosity 021. Our final sample for fitting A/bh — i contains 
36 black hole masses. 

Our fits to A/bh — o and A/bh — i assume a single-index power law 
as the functional form of both relations, following the convention most 
commonly used in prior studies. Specifically, we define the A/bh — o" 
relationship to be logj^Q(MBH/ A/q) = a ^r /3 logjQ(cr/200 km s^^), 
and the A/bh — L relationship to be logjo(A/BH/ A/©) = a + 
/9 logj^Q(Lv/10^^LQ). For each relationship, we follow the method of 



Tremaine et al. (2002)1^10 fit for a and /3. We minimize the quantity 

2 _ ■A (Mbh,. -a - Pojf ,. 

X - 2 , 2 , a2 2 

where is the measurement error in o, eu is the measurement error in 
Mbh, and eo is the intrinsic scatter in the A/bh — relation. We set eq 
such that per degree of freedom is unity after minimization. The 68% 
confidence intervals for a and /3 correspond to the maximum range of a 
and P for which - xLin < 1- 

We list our best-fit values of a, /3, and eo for various sub-samples in 
Table m We list our full sample of galaxies in Table [4] Graham et al. 
(20 11)^ recently compiled a sample of 64 galaxies with dynamical black 
hole measurements and found a power-law index of 5.13 for the A/bh — c 
relation, very similar to the index we report herein. However, only 52 
galaxies appear in their sample as well as ours. Of the 13 galaxies that 
appear only i n our sample, seven use precise maser-based measurements 
of Mbh'^^^, and three are BCGs. Additionally, 1 1 galaxies have received 
updated measurements of A/bh since the compilation of Graham et al. 
(2011), using stellar orbit models with dark matter halos and larger orbit 
libraries'^. 
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Table 1 | Spectroscopic observations 



Galaxy 


Instrument 


UTDate 


ATexp 




A 


AA 




Ar 


PSF FWHM 










(s) 


(nm) 


(nm) 


(arcsec) 


(arcsec) 


(arc sec) 


NGC 3842 


GMOS 


April 27, 2003 


5 


1200 


755-949 


0.25 


3.8 


0.2 


0.4 


NGC 3842 


OSIRIS 


May 8-10, 2010 


16 


900 


1473-1803 


0.72 


0.7 


0.05 


0.08 


NGC 3842 


VIRUS-P 


March 7-8, 2011 


3 


1200 


358-589 


0.52 


35.3 


4.1 


2.0 


NGC 4889 


GMOS 


March 13, 2003 


5 


1200 


755-949 


0.25 


3.8 


0.2 


0.4 



Notes: A^exp is the number of science exposures recorded, fexp is the integration time per exposure. AA is the median instrumental resolution (FWHM). 
r-max is the maximum radius of usable data, with respect to the center of the galaxy. Ar is the instrumental spatial scale, equal to the angular diameter 
of one lenslet or fiber. 



Table 2 | IVIodeis of NGC 4889 



Quadrant 


Mbh 


A/BH,min 


A^BH,max 












(Me) 


(Mq) 




(Mq/Lq) 


(Mo/Lq) 


northeast 


1.7 X 10^° 


1.0 X 10^° 


2.5 X 10^° 


6.1 


4.6 


7.3 


southeast 


2.6 X 10^" 


2.0 X 10^° 


3.2 X 10^° 


5.6 


4.2 


6.7 


northwest 


2.7 X 10^° 


1.6 X 10^° 


3.7 X 10^° 


5.8 


4.4 


7.0 


southwest 


9.8 X 10^ 


5.5 X 10^ 


1.7 X 10^° 


6.6 


5.3 


7.6 


east 


2.9 X 10^" 


2.1 X 10^" 


3.4 X 10"' 


5.4 


4.5 


6.4 


west 


1.2 X 10^° 


6.5 X 10^ 


2.0 X 10^° 


6.4 


5.2 


7.4 


recentered 


1.5 X 10^° 


8.7 X 10^ 


2.4 X 10^° 


6.5 


5.4 


7.3 



Notes: The "east" and "west" trials used LOSVDs from spectra that were binned symmetrically over the north and south sides of the galaxy. The "re- 
centered" trial added an artificial position offset to the LOSVDs, such that the maximum velocity dispersion was placed at the center of the gravitational 
potential. MBH,min, Mbh, max, M,/I/, min, and Mt/L,inax represent 68% confidence limits. M*/L corresponds to il-band photometry. 



Table 3 | Fits to Mbh - o- and Mbh - L 



Relationship 


Sample 


A^gal 


a 




eo 


Mbh - 0- 


all 


65 


8.29 ±0.06 


5.12 ±0.36 


0.43 


Mbh - a 


early-type 


45 


8.38 ±0.06 


4.53 ±0.40 


0.38 


Mbh - a- 


late-type 


20 


7.97 ±0.22 


4.58 ± 1.25 


0.44 


Mbh - a- 


ML 


36 


8.43 ± 0.07 


4.66 ± 0.43 


0.38 


Mbh - L 


ML 


36 


9.16 ±0.11 


1.16 ±0.14 


0.50 



Notes: TVgai is the number of galaxies in the sample. ML refers to the sample of early-type galaxies with reliable spheroid luminosity measurements. 
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Table 4 | Galaxies with dynamical measurements of Mbh 



Galaxy Type D a logio(I/v,sph) Mbh AfeHimin AfBH,max Method Ref. 



(Mpc) (kms-i) logio(L0) (Mq) (Mq) (Mq) 



Milky Way " 


Sbc 


0.008 


103 ± 20 




4.1 X 


10*^ 


3.5 


X 


10® 


4.7 X 


10® 


stars 


|9|41|42| 


A1836-BCG 


E 


157.5 


288 ± 14 


11.26 


3.9 X 


10« 


3.3 


X 


10^ 


4.3 X 


lo'-* 


gas 


IWl 


A3565-BCG '' 


E 


54.4 


322 ± 16 


11.24 


1.4 X 


10» 


1.2 


X 


10^ 


1.7 X 


10^ 


gas 


m 


Circinus 


Sb 


4.0 


158 ± 18 




1.7 X 


10^ 


1.4 


X 


10® 


2.1 X 


10® 


masers 


IWl 


IC1459 " 


E4 


30.9 


315 ± 16 


10.96 


2.8 


X 


10« 


1.6 


X 


10» 


3.9 X 


10^ 


stars 


EEl 


N221 (M32) 


E2 


0.86 


75 ±3 


8.66 


3.1 X 


10'' 


2.5 


X 


10® 


3.7 X 


10® 


stars 


EE3 


N224 (M31) 


Sb 


0.80 


160 ±8 




1.5 X 


10** 


1.2 


X 


10** 


2.4 X 


10* 


stars 


EEH 


N524 


SO 


23.3 


235 ± 12 


10.67 


8.3 X 


10« 


7.9 


X 


10** 


9.2 X 


10** 


stars 


M 


N821 


E4 


25.5 


209 ± 10 


10.43 


1.8 


X 


10** 


1.0 


X 


10** 


2.6 X 


10** 


stars 


m 


N1023 


SBO 


12.1 


205 ± 10 


10.18 


14.6 X 


10^ 


4.1 


X 


10^ 


5.1 X 


10^ 


stars 


lag 


N1068 (M77) 


SB 


15.4 


151 ±7 




8.6 X 


10^ 


8.3 


X 


10® 


8.9 X 


10® 


masers 


011 


N1194 ''•'^ 


SAO+ 


55.5 


U8tll 




6.8 


X 


10^ 


6.5 


X 


10^ 


7.1 X 


10^ 


masers 


m 


N1300 


SB(rs)bc 


20.1 


218 ± 10 




7.1 X 


10^ 


3.6 


X 


10^ 


1.4 X 


10** 


gas 


Em 


N1316 


E 


18.6 


226 ± 11 


11.26 


1.5 X 


10* 


1.24 X 


10** 


1.75 X 


10* 


stars 


m 


N1332 


SO 


22.3 


328 ± 16 


10.14 


1.45 X 


10^ 


1.25 X 


10^ 


1.65 X 


10^ 


stars 


m 


N1399 "-^ 


El 


21.1 


296 ± 15 


10.78 


5.1 X 


10** 


4.4 X 


10** 


5.8 X 


10** 


stars 


Em 


N1399 "'f' 


El 


21.1 


296 ± 15 


10.78 


1.3 X 


10« 


6.4 X 


10** 


1.8 X 


10^ 


stars 


Em 


N2273 ''•'^ 


SB(r)a 


26.8 


14411^ 




7.8 X 


10^ 


7.4 X 


10® 


8.2 X 


10® 


masers 


m 


N2549 


SO 


12.3 


145 ±7 


9.60 


1.4 X 


10^ 


1.0 X 


10^ 


1.47 X 


10^ 


stars 


EH 


N2748 


Sc 


24.9 


115 ±5 




4.7 X 


10^ 


8.6 X 


10® 


8.5 X 


10^ 


gas 


Em 


N2787 


SBO 


7.9 


189 ±9 




4.3 X 


10^ 


3.8 X 


10^ 


4.7 X 


10^ 


gas 


Em 


N2960 ''•'^ 


Sa 


75.3 






1.21 X 


10^ 


1.16 X 


10^ 


1.26 X 


10^ 


masers 


m 


N3031 (M81) 


Sb 


4.1 


143 ± 7 




8.0 


X 


10^ 


6.9 


X 


10^ 


1.0 X 


10** 


gas 


Em 


N3115 


SO 


10.2 


230 ± 11 


10.40 


9.6 


X 


10** 


6.7 


X 


10** 


1.5 X 


10^ 


stars 


Em 


N3227 


SBa 


17.0 


133 ± 12 




1.5 


X 


10^ 


7.0 


X 


10® 


2.0 X 


10^ 


stars 


Em 


N3245 


SO 


22.1 


205 ± 10 




2.2 


X 


10** 


1.7 


X 


10** 


2.7 X 


10** 


gas 


Em 


N3368 ^ 


SAB(rs)ab 


10.4 


122tll 




7.5 


X 


10® 


6.0 


X 


10® 


9.0 X 


10® 


stars 


m 


N3377 


E6 


11.7 


145 ±7 


9.98 


1.9 


X 


10** 


9.0 


X 


10^ 


2.9 X 


10** 


stars 


m 


N3379 (M105) 


EO 


11.7 


206 ± 10 


10.37 


4.6 


X 


10** 


3.4 


X 


10** 


5.7 X 


10** 


stars 


m 


N3384 


El 


11.7 


143 ± 7 


9.90 


1.1 


X 


10^ 


6.0 


X 


10® 


1.6 X 


10^ 


stars 


m 


N3393 '' 


SB(rs) 


53.6 


148 ± 10 




3.3 


X 


10^ 


3.1 


X 


10^ 


3.5 X 


10^ 


masers 


M 


N3489 s 


SAB(rs)0+ 


12.1 


loolli 




6.0 


X 


10*^ 


5.2 


X 


10® 


6.8 X 


10® 


stars 


M 


N3585 


SO 


21.2 


213 ± 10 


10.65 


3.4 


X 


10** 


2.8 


X 


10** 


4.9 X 


10** 


stars 


Em 


N3607 


El 


19.9 


229 ± 11 




1.2 


X 


10** 


7.9 


X 


10^ 


1.6 X 


10** 


stars 


Em 


N3608 


El 


23.0 


182 ±9 


10.35 


4.7 


X 


10** 


3.7 


X 


10** 


5.7 X 


10** 


stars 


m 


N3842 ' 


E 


98.4 


270 ± 14 


11.20 


9.7 


X 


10« 


7.2 


X 


10^ 


1.27 X 


10^° 


stars 


Em 


N3998 


SO 


14.9 


305 ± 15 




2.4 


X 


10** 


6.2 


X 


10^ 


4.5 X 


10** 


gas 


N4026 


SO 


15.6 


180 ±9 


9.86 


2.1 


X 


10** 


1.7 


X 


10** 


2.8 X 


10** 


stars 


Em 


N4258 


SABbc 


7.2 


115 ± 10 




3.78 X 


10^ 


3.77 X 


10^ 


3.79 X 


10^ 


masers 


Em 


N4261 


E2 


33.4 


315 ± 15 


11.02 


5.5 


X 


10** 


4.3 


X 


10** 


6.6 X 


10** 


gas 


Em 


N4291 


E2 


25.0 


242 ± 12 


10.20 


9.2 


X 


10** 


6.3 


X 


10** 


1.21 X 


10^ 


stars 


m 


N4342 


SO 


18.0 


225 ± 11 




3.6 


X 


10** 


2.4 


X 


10** 


5.6 X 


10** 


stars 


Em 


N4374 (M84) 


El 


17.0 


296 ± 14 


10.91 


8.5 


X 


10** 


7.7 


X 


10** 


9.4 X 


10** 


gas 


m 


N4388 ''•'^ 


SA(s)b 


19.8 


107t? 




8.8 


X 


10® 


8.6 


X 


10® 


9.0 X 


10® 


masers 


m 


N4459 


E2 


17.0 


167 ±8 


10.36 


7.4 


X 


10^ 


6.0 


X 


10^ 


8.8 X 


10^ 


gas 


Em 


N4473 


E4 


17.0 


190 ±9 


10.39 


1.0 


X 


10« 


5.0 


X 


10^ 


1.5 X 


10** 


stars 


m 


N4486 (M87) ''■^ 


El 


17.0 


3241?^ 


11.10 


6.3 


X 


10» 


5.9 


X 


10^ 


6.6 X 


10^ 


stars 


El 


N4486A 


E2 


17.0 


Ill ± 5 


9.41 


1.3 


X 


10^ 


9.0 


X 


10® 


1.8 X 


10^ 


stars 


Em 


N4564 


SO 


17.0 


162 ±8 




9.4 


X 


10^ 


6.8 


X 


10^ 


1.2 X 


10** 


stars 


m 


N4594 (Ml 04) 


Sa 


10.3 


240 ± 12 




5.3 


X 


10** 


4.74 X 


10** 


6.08 X 


10** 


stars 


m 



Continued on next page. 
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Table 4, continued 



Galaxy Type D a logio(I/v,aph) Mbh MBH,min MBH,max Method Ref. 







(Mpc) 


(kms-^) 


logio(L0) 


(Me) 




(Mq) 






N4596 


SBO 


18.0 


136 ±6 




8.4 X 10^ 


5.9 X 10^ 


1.2 X 10* 


gas 


EEl 


N4649 (M60) 


E2 


16.5 


341 ± 17 


10.99 


4.7 X 10^ 


3.7 X 10^ 


5.8 X lO'* 


stars 


w\ 


N4697 


E6 


12.4 


177 ±8 


10.45 


2.0 X 10** 


1.8 X 10* 


2.2 X 10* 


stars 


m 


N4736 (M94) 


Sab 


4.9 


112 ±6 




6.68 X 10^ 


5.14 X 10*^ 


8.22 X lO'^ 


stars 


m 


N4826 (M64) 


Sab 


6.4 


96 ±5 




1.36 X 10^ 


1.02 X 10® 


1.71 X lO'^ 


stars 


m 


N4889 


E 


103.2 


347 ± 17 


11.42 


2.1 X 10^° 


5.5 X 10^ 


3.7 X 10^" 


stars 


Em 


N5077 


E3 


44.9 


222 ± 11 


10.75 


8.0 X 10* 


4.7 X 10* 


1.3 X lO'' 


gas 


N5128 " 


SO/E 


4.4 


150 ±7 


10.66 


3.0 X 10** 


2.8 X 10* 


3.4 X 10* 


stars 


EE!] 


N5128 = 


SO/E 


4.4 


150 ±7 


10.66 


7.0 X 10^ 


3.2 X 10^ 


8.3 X 10^ 


stars 


13221 


N5576 


E3 


27.1 


183 ±9 


10.44 


1.8 X 10* 


1.4 X 10* 


2.1 X 10* 


stars 


M^ii 


N5845 


E3 


28.7 


234 ± 11 


9.84 


5.4 X 10* 


3.7 X 10* 


7.1 X 10* 


stars 


[U 


N6086 


E 


139.1 


318 ± 16 


11.18 


3.8 X 10^ 


2.6 X 10^ 


5.5 X 10^ 


stars 


EH 


N6251 


El 


106.0 


290 ± 14 




6.0 X 10* 


4.0 X 10* 


8.0 X 10* 


gas 


EED 


N6264 ''''' 


S 


145.4 


158tli 




3.03 X 10^ 


2.99 X 10^ 


3.08 X 10^ 


masers 


[U 


N6323 


Sab 


110.5 


15811 




9.8 X lO'' 


9.7 X lO'' 


9.9 X lO'* 


masers 


[m 


N7052 


E3 


70.9 


266 ± 13 


10.87 


4.0 X 10* 


2.4 X 10* 


6.8 X 10* 


gas 




N7457 


SO 


14.0 


67 ±3 


9.42 


1.0 X 10^ 


4.0 X lO'' 


1.6 X 10^ 


stars 


mi 


N7582 


SBab 


22.3 


156 ± 19 




5.5 X 10^ 


4.4 X 10^ 


7.1 X 10^ 


gas 


|9l73l 


U3789 


SA(r)ab 


48.4 


107111 




1.08 X 10^ 


1.03 X 10^ 


1.14 X 10'' 


masers 





Notes: D is the distance to the galaxy, assuming Hq — 70kms Mpc . Lv.sph is the V-band stellar luminosity of the galaxy's spheroidal 
component. A/BH,miii and A/BH,max are the upper and lower bounds of the 68% confidence interval in A/bh. 

" The literature contains a large number of estimates for the velocity dispersion of our Galaxy's bulge, using different kinematic tracers at different 
radii. We use the radially averaged measurement of cr = 103 ± 20 km s^^ from Tremaine et al. (2002) 

G09 use the convention Ho = 70kms^'^ Mpc^^ for all galaxies whose distances are derived from systemic velocity measurements. To match this 
convention, we have adjusted the distance to several galaxies in our updated sample. All distance adjustments yield corresponding adjustments to 
A^bh, A^BHimin, and A/bh, max, such that A/bh oc D. Our reported 68% confidence intervals for A/bh do not include uncertainties in the distance 
measurements. 

We derive an effective velocity dispersion of 315 kms^^ for IC 1459, using spatially resolved major-axis kinematics from Cappellari et al. (2002)1^, 
at radii of 0.8 - 25.1 arcseconds. For r < 0.8 arcseconds, stellar motions are directly influenced by the central black hole. 

Maser-based black hole masses for several galaxies are presented in Greene et al. (2010) and Kuo et al. (201 1) We use the velocity dispersions 
presented in Greene et al. (2010). For consistency with the rest of our sample, we use the black hole masses from Kuo et al. (2011), which agree 
with the values in Greene et al. (2010) but do not include distance uncertainties in the overall uncertainty for A/bh. Braatz et al. (20 10) '^provide an 
updated distance and black hole mass for UGC 3789, which are consistent with the values we adopt from Kuo et al. (201 1). 

Following G09, our sample includes two distinct measurements for NGC 1399 and also for NGC 5128. We weight each of these measurements by 
50% when performing fits to MBti{o') and A^BH(i). In Figurejsjwe only plot the more recent measurement for each galaxy. 

^ We derive an effective velocity dispersion of 296kms~^ for NGC 1399, using spatially resolved measurements from Graham et al. (1998)1^ and 
Gebhardt et al. (2007)'^ at radii of 0.6 - 41 arcseconds. For r < 0.6 arcseconds, stellar motions are directly influenced by the central black hole. 
® NGC 3368 and NGC 3489 were not included in our sample until late in the editorial process, and so do not contribute toward the fits reported in the 
main letter or in Table 3. Still, we view the reported black hole masses for NGC 3368 and NGC 3489 as valid measurements that should be included 
in future studies. For the velocity dispersions, we use the average values and outer limits of the various measurements reported by Nowak et al. (2010) 
EH. Including NGC 3368 and NGC 3489, the A/bh - cr parameters (a, /3, eo) are (8.28 ± 0.06, 5.13 ± 0.34, 0.42) for 67 galaxies, (8.38 ± 0.06, 
4.57 ± 0.36, 0.36) for 46 early-type galaxies, and (8.00 ± 0.21, 4.76 ± 1.15, 0.43) for 21 late-type galaxies. These parameters are consistent with our 
reported values for the sample of 65 galaxies. 

The literature contains two inconsistent estimates of the stellar luminosity in NGC 3607: My = —21.62 in G09, and My ~ —19.88 in Lauer et al. 
(2007)G3. 

' We derive an effective velocity dispersion of 270 km s^^ for NGC 3842, using measurements from GMOS and VIRUS-P at radii of 1.2 - 29.8 
arcseconds. For r < 1.2 arcseconds, stellar motions are directly influenced by the central black hole. 

For M87, we use the updated velocity dispersion of 324 km s^^ from Gebhardt et al. (201 1) based on measurements at radii of 2.1 - 100 arcseconds. 
For r < 2.1 arcseconds, stellar motions are directly influenced by the central black hole. 

The literature contains two inconsistent estimates of the bulge stellar luminosity in NGC 4564: My = —19.60 in G09, and My = —20.26 in Lauer 
et al. (2007)P 

' We derive an effective velocity dispersion of 341 kms^^ for M60, using spatially resolved measurements from Pinkney et al. (2003)'^ at radii of 2.2 
- 44 arcseconds. For r < 2.2 arcseconds, stellar motions are directly influenced by the central black hole. 

™ Our quoted value of A/bh for NGC 4889 is the median of the 68% confidence interval 0.7 — 3.4 x 10^° A/q. We derive an effective velocity 
dispersion of 347kms~^, using kinematics from Loubser et al. (2008) at radii of 2.0 - 22.9 arcseconds. The inner cutoff of 2.0 arcseconds excludes 
the asymmetric velocity dispersion peak and possible stellar torus. 
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